Integrating dynamical modeling and phylogeographic inference to characterize global influenza circulation.

Global seasonal influenza circulation involves a complex interplay between local (seasonality, demography, host immunity) and global factors (international mobility) shaping recurrent epidemic patterns. No studies so far have reconciled the two spatial levels, evaluating the coupling between national epidemics, considering heterogeneous coverage of epidemiological and virological data, integrating different data sources. We propose a novel combined approach based on a dynamical model of global influenza spread (GLEAM), integrating high-resolution demographic and mobility data, and a generalized linear model of phylogeographic diffusion that accounts for time-varying migration rates. Seasonal migration fluxes across global macro-regions simulated with GLEAM are tested as phylogeographic predictors to provide model validation and calibration based on genetic data. Seasonal fluxes obtained with a specific transmissibility peak time and recurrent travel outperformed the raw air-transportation predictor, previously considered as optimal indicator of global influenza migration. Influenza A subtypes supported autumn-winter reproductive number as high as 2.25 and an average immunity duration of 2 years. Similar dynamics were preferred by influenza B lineages, with a lower autumn-winter reproductive number. Comparing simulated epidemic profiles against FluNet data offered comparatively limited resolution power. The multiscale approach enables model selection yielding a novel computational framework for describing global influenza dynamics at different scales - local transmission and national epidemics vs. international coupling through mobility and imported cases. Our findings have important implications to improve preparedness against seasonal influenza epidemics. The approach can be generalized to other epidemic contexts, such as emerging disease outbreaks to improve the flexibility and predictive power of modeling.

Phylogenetics and Population Genetics of the Asian House Shrew, Suncus murinus-S. montanus Species Complex, Inferred From Whole-Genome and Mitochondrial DNA Sequences, with Special Reference to the Ryukyu Archipelago, Japan.

The house shrew (Suncus murinus-S. montanus species complex) colonized regions across southern Asia and the Indian Ocean following human activity. The house shrew is distributed on islands of the Ryukyu Archipelago, the southernmost part of Japan, but the evolutionary history of the shrew on those islands and possible associations between these populations and humans remain unknown. In this study, we conducted phylogenetic and population genetic analyses based on both nuclear and mitochondrial genome sequences of house shrews. Phylogenetic analyses based on mitochondrial cytochrome b (cytb) sequences revealed that shrews from the Ryukyu Archipelago showed strong genetic affinity to Vietnamese and southern Chinese shrews. Demographic analyses of cytb sequences indicated a rapid population expansion event affecting the haplotype group in Vietnam, southern China, and the Ryukyu Archipelago 3300-7900 years ago. Furthermore, gene flow between Ryukyu (Yonaguni Island) and Taiwan and between Ryukyu and Vietnam inferred from f4 statistics of the nuclear genomes suggested repeated immigration to Ryukyu in recent years. The present study demonstrates that the Nagasaki population has a different origin from the Ryukyu population. These findings elucidate the complex pattern of genetic admixture in house shrews and provide insights into their evolutionary history.

Revision of the 'Acanthephyra purpurea' species complex (Crustacea: Decapoda), with an emphasis on species diversification in the Atlantic.

We inventoried all nine species of the 'Acanthephyra purpurea' complex, one of the most abundant and cosmopolitan group of mesopelagic shrimps. We used 119 specimens at hand and genetic data for 124 specimens from GenBank and BOLD. Phylogenetic analysis of four genes (COI, 16S, NaK, and enolase) showed that the 'Acanthephyra purpurea' complex is polyphyletic and encompasses two species groups, 'A. purpurea' (mostly Atlantic) and 'A. smithi' (Indo-West Pacific). The 'A. purpurea' species group consists of two major molecular clades A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa. Molecular data suggest that hitherto accepted species A. acanthitelsonis, A. pelagica, and A. sica should be considered as synonyms. The Atlantic is inhabited by at least two cryptic genetic lineages of A. pelagica and A. quadrispinosa. Morphological analyses of qualitative and quantitative (900 measurements) characters resulted in a tabular key to species and in a finding of four evolutionary traits. Atlantic species showed various scenarios of diversification visible on mitochondrial gene level, nuclear gene level, and morphological level. We recorded and discussed similar phylogeographic trends in diversification and in distribution of genetic lineages within two different clades: A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa.

A cryptic radiation of Caribbean sea slugs revealed by integrative analysis: Cyerce 'antillensis' (Sacoglossa: Caliphyllidae) is six distinct species.

Integrative studies have revealed cryptic radiations in several Caribbean lineages of heterobranch sea slugs, raising questions about the evolutionary mechanisms that promote speciation within the tropical Western Atlantic. Cyerce Bergh, 1871 is a genus comprising 12 named species in the family Caliphyllidae that lack the photosynthetic ability of other sacoglossans but are noted for vibrant colours on the large cerata (dorsal leaf-like appendages) that characterize many species. Two species are widely reported from the Caribbean: Cyerce cristallina (Trinchese, 1881) and Cyerce antillensis Engel, 1927. Here, we present an integrative assessment of diversity in Caribbean Cyerce. Four methods of molecular species delimitation supported seven species in samples from the Caribbean and adjacent subtropical Western Atlantic. Six delimited species formed a monophyletic lineage in phylogenetic analyses but were > 9% divergent at the barcoding COI locus and could be differentiated using ecological, reproductive and/or morphological traits. We redescribe C. antillensis, a senior synonym for the poorly known Cyerce habanensis Ortea & Templado, 1988, and describe five new species. Evolutionary shifts in algal host use, penial armature and larval life history might have acted synergistically to promote the rapid divergence of endemic species with restricted distributions in this radiation, substantially increasing global diversity of the genus.

Population Variation and Phylogeography of Cherry Blossom (Prunus conradinae) in China.

Prunus conradinae (subgenus Cerasus, Rosaceae) is a significant germplasm resource of wild cherry blossom in China. To ensure the comprehensiveness of this study, we used a large sample size (12 populations comprising 244 individuals) which involved the fresh leaves of P. conradinae in Eastern, Central, and Southwestern China. We combined morphological and molecular evidence (three chloroplast DNA (cpDNA) sequences and one nuclear DNA (nr DNA) sequence) to examine the population of P. conradinae variation and differentiation. Our results revealed that Central, East, and Southwest China are important regions for the conservation of P. conradinae to ensure adequate germplasm resources in the future. We also found support for a new variant, P. conradinae var. rubrum. We observed high genetic diversity within P. conradinae (haplotype diversity [Hd] = 0.830; ribotype diversity [Rd] = 0.798), with novel genetic variation and a distinct genealogical structure among populations. There was genetic variation among populations and phylogeographic structure among populations and three geographical groups (Central, East, and Southwest China). The genetic differentiation coefficient was the lowest in the Southwest region and the gene exchange was obvious, while the differentiation was obvious in Central China. In the three geographic groups, we identified two distinct lineages: an East China lineage (Central China and East China) and a Southwest China lineage ((Central China and Southwest China) and East China). These two lineages originated approximately 4.38 million years ago (Mya) in the early Pliocene due to geographic isolation. P. conradinae expanded from Central China to East China at 3.32 Mya (95% HPD: 1.12-5.17 Mya) in the Pliocene. The population of P. conradinae spread from East China to Southwest China, and the differentiation time was 2.17 Mya (95% (HPD: 0.47-4.54 Mya), suggesting that the population of P. conradinae differentiated first in Central and East China. The population of P. conradinae experienced differentiation from Central China to Southwest China around 1.10 Mya (95% HPD: 0.11-2.85 Mya) during the early Pleistocene of the Quaternary period. The southeastern region of East China, near Mount Wuyi, likely serves as a refuge for P. conradinae. This study establishes a theoretical foundation for the classification, identification, conservation, and exploitation of germplasm resources of P. conradinae.

Genetic differentiation and phylogeography of Erythroneurini (Hemiptera, Cicadellidae, Typhlocybinae) in the southwestern karst area of China.

Erythroneurini is the largest tribe of the microleafhopper subfamily Typhlocybinae. Most prior research on this tribe has focused on traditional classification, phylogeny, and control of agricultural pests, and the phylogeography of the group remains poorly understood. In this study, the mitochondrial genomes of 10 erythroneurine species were sequenced, and sequences of four genes were obtained for 12 geographical populations of Seriana bacilla. The new sequence data were combined with previously available mitochondrial DNA sequence data and analyzed using Bayesian and Maximum-Likelihood-based phylogenetic methods to elucidate relationships among genera and species and estimate divergence times. Seriana was shown to be derived from within Empoascanara. Phylogeographic and population genetic analysis of the endemic Chinese species Seriana bacilla suggest that the species diverged about 54.85 Mya (95% HPD: 20.76-66.23 million years) in the Paleogene period and that population divergence occurred within the last 14 million years. Ancestral area reconstruction indicates that Seriana bacilla may have originated in the central region of Guizhou, and geographical barriers are the main factors affecting gene flow among populations. Ecological niche modeling using the MaxEnt model suggests that the distribution of the species was more restricted in the past but is likely to expand in the future years 2050 and 2070.

DNA barcoding, phylogeography and evolutionary dynamics of Chrysichthys auratus.

This study embarked on an exploration into the genetic structure and evolutionary history of the Chrysichthys auratus species, leveraging PCR amplification, phylogenetic trees, and haplotype networks. Specific DNA segments were successfully amplified and visualized through electrophoresis. Newly obtained sequences were Bank into GenBank and given accession numbers (OR730807-OR730808-OR730809). The Neighbor-Joining method provided insights into the evolutionary relationships among taxa, further augmented by bootstrap values and the Tamura 3-parameter method. A comprehensive geographical haplotype network showcased pronounced genetic differentiation, especially between remote populations. Nonetheless, shared haplotypes between proximate regions indicated either ancestral genetic connections or ongoing gene flow. Employing the COI-DNA barcodes, an in-depth understanding of intra- and inter-populational genetic diversity was achieved. The study's findings unravel the intricate genetic landscape and evolutionary dynamics of C. auratus, offering novel perspectives into its demographic history across its vast native habitat.

Phylodynamic analysis of foot-and-mouth disease virus evolution in Mar Chiquita, Argentina.

Foot-and-mouth disease is a highly contagious disease affecting cloven-hoofed animals, resulting in considerable economic losses. Its causal agent is foot-and-mouth disease virus (FMDV), a picornavirus. Due to its error-prone replication and rapid evolution, the transmission and evolutionary dynamics of FMDV can be studied using genomic epidemiological approaches. To analyze FMDV evolution and identify possible transmission routes in an Argentinean region, field samples that tested positive for FMDV by PCR were obtained from 21 farms located in the Mar Chiquita district. Whole FMDV genome sequences were obtained by PCR amplification in seven fragments and sequencing using the Sanger technique. The genome sequences obtained from these samples were then analyzed using phylogenetic, phylogeographic, and evolutionary approaches. Three local transmission clusters were detected among the sampled viruses. The dataset was analyzed using Bayesian phylodynamic methods with appropriate coalescent and relaxed molecular clock models. The estimated mean viral evolutionary rate was 1.17 × 10- 2 substitutions/site/year. No significant differences in the rate of viral evolution were observed between farms with vaccinated animals and those with unvaccinated animals. The most recent common ancestor of the sampled sequences was dated to approximately one month before the first reported case in the outbreak. Virus transmission started in the south of the district and later dispersed to the west, and finally arrived in the east. Different transmission routes among the studied herds, such as non-replicating vectors and close contact contagion (i.e., aerosols), may be responsible for viral spread.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics allowing to unveil the dispersal dynamics and evaluate the capacity to spread among hosts. Heterogeneous sampling intensity of genomic sequences can however impact the accuracy of dispersal insights gained through phylogeographic inference. In our study, we implement a simulation framework to evaluate the robustness of three dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance signal metric - to the sampling effort. Our results reveal that both the diffusion coefficient and isolation-by-distance signal metrics appear to be robust to the number of samples considered for the phylogeographic reconstruction. We then use these two dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of isolation-by-distance patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Conserved DNA sequence analysis reveals the phylogeography and evolutionary events of Akebia trifoliata in the region across the eastern edge of the Tibetan Plateau and subtropical China.

BACKGROUND: The eastern edge of the Qinghai‒Tibet Plateau (QTP) and subtropical China have various regions where plant species originate and thrive, but these regions have been the focus of very few integrative studies. Here, we elucidated the phylogeographic structure of a continuous and widespread Akebia trifoliata population across these two regions. RESULTS: Sixty-one populations consisting of 391 genotypes were examined to assess population diversity and structure via network distribution analysis, maximum likelihood phylogenetic tree reconstruction, divergence time estimation, demographic history inference, and ancestral area reconstruction of both conserved internal transcribed spacer (ITS) and chloroplast (rps16) DNA sequences. The results showed that the ITS region was more variable than the rps16 region and could be suitable for studying intraspecific phylogeography. The A. trifoliata population displayed high genetic diversity, genetic differentiation and obvious phylogeographical structure, possibly originating on the eastern QTP, expanding during the last glacial-interglacial cycle, diverging in the early Pleistocene and middle Pleistocene, and extensively migrating thereafter. The migration route from west to east along rivers could be largely responsible for the long-distance dispersal of this species, while three main refuges (Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) with multiple ice shelters facilitated its wide distribution. CONCLUSIONS: Our results suggested that the from west to east long migration accompanying with the minor short reciprocal migration in the south-north direction, and the three main refuges (the Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) contributed to the extant geographical distribution of A. trifoliata. In addition, this finding also strongly reduced the discrepancy between glacial contraction and postglacial expansion and the in situ survival hypothesis by simultaneously considering the existence of many similar climate-related ecological niches and migration influences.

Environmental DNA unveils deep phylogeographic structure of a freshwater fish.

Phylogeography bears an important part in ecology and evolution. However, current phylogeographic studies are largely constrained by limited numbers of individual samples. Using an environmental DNA (eDNA) assay for phylogeographic analyses, this study provides detailed information regarding the history of Siberian stone loach Barbatula toni, a primary freshwater fish across the whole range of Hokkaido, Japan. Based on an eDNA metabarcoding on 293 river water samples, we detected eDNA from B. toni in 189 rivers. A total of 51 samples, representing the entire island, were then selected from the B. toni eDNA-positive sample set for the subsequent analyses. To elucidate the phylogeographic structure of B. toni, newly developed eDNA metabarcoding primers (Barba-cytb-F/R) were applied to these samples, specifically targeting their haplotypic variation in cytochrome b. After a bioinformatic processing to mitigate haplotypic false positives, a total of 50 eDNA haplotypes were identified. Two regionally restricted, genetically distinct lineages of the species were revealed as a result of phylogeographic analyses on the haplotypes and tissue-derived DNA from B. toni. According to a molecular clock analysis, they have been genetically isolated for at least 1.5 million years, suggesting their ancient origin and colonisation of Hokkaido, presumably in the glacial periods. These results demonstrate how freshwater fishes can alter their distributions over evolutionary timescales and how eDNA assay can deepen our understanding of phylogeography.

Survival in nunatak and peripheral glacial refugia of three alpine plant species is partly predicted by altitudinal segregation.

Mountain biota survived the Quaternary cold stages most probably in peripheral refugia and/or ice-free peaks within ice-sheets (nunataks). While survival in peripheral refugia has been broadly demonstrated, evidence for nunatak refugia is still scarce. We generated RADseq data from three mountain plant species occurring at different elevations in the southeastern European Alps to investigate the role of different glacial refugia during the Last Glacial Maximum (LGM). We tested the following hypotheses. (i) The deep Piave Valley forms the deepest genetic split in the species distributed across it, delimiting two peripheral refugia. (ii) The montane to alpine species Campanula morettiana and Primula tyrolensis survived the LGM in peripheral refugia, while high-alpine to subnival Saxifraga facchinii likely survived in several nunatak refugia. (iii) The lower elevation species suffered a strong population decline during the LGM. By contrast, the higher elevation species shows long-term stability of population sizes due to survival on permanently ice-free peaks and small population sizes at present. We found peripheral refugia on both sides of the Piave Valley, which acted as a major genetic barrier. Demographic modelling confirmed nunatak survival not only for S. facchinii but also for montane to alpine C. morettiana. Altitudinal segregation influenced the species' demographic fluctuations, with the lower elevation species showing a significant population increase at the end of the LGM, and the higher elevation species either showing decrease towards the present or stable population sizes with a short bottleneck. Our results highlight the role of nunatak survival and species ecology in the demographic history of mountain species.

The Andes as a semi-permeable geographical barrier: Genetic connectivity between structured populations in a widespread spider.

Geographical barriers like mountain ranges impede genetic exchange among populations, promoting diversification. The effectiveness of these barriers in limiting gene flow varies between lineages due to each species' dispersal modes and capacities. Our understanding of how the Andes orogeny contributes to species diversification comes from well-studied vertebrates and a few arthropods and plants, neglecting organisms unable to fly or walk long distances. Some arachnids, such as Gasteracantha cancriformis, have been hypothesized to disperse long distances via ballooning (i.e. using their silk to interact with the wind). Yet, we do not know how the environment and geography shape its genetic diversity. Therefore, we tested whether the Andes contributed to the diversification of G. cancriformis acting as an absolute or semi-permeable barrier to genetic connectivity between populations of this spider at opposite sides of the mountain range. We sampled thousands of loci across the distribution of the species and implemented population genetics, phylogenetic, and landscape genetic analyses. We identified two genetically distinct groups structured by the Central Andes, and a third less structured group in the Northern Andes that shares ancestry with the previous two. This structure is largely explained by the altitude along the Andes, which decreases in some regions, possibly facilitating cross-Andean dispersal and gene flow. Our findings support that altitude in the Andes plays a major role in structuring populations in South America, but the strength of this barrier can be overcome by organisms with long-distance dispersal modes together with altitudinal depressions.

Las barreras geográficas como las cordilleras montañosas impiden el intercambio genético entre poblaciones, promoviendo la diversificación. La efectividad de estas barreras para limitar el flujo genético varía entre linajes debido a los modos y capacidades de dispersión de cada especie. Nuestra comprensión de cómo la orogenia de los Andes contribuye a la diversificación de especies proviene de vertebrados y algunos artrópodos y plantas bien estudiados, descuidando a los organismos incapaces de volar o caminar grandes distancias. Se ha hipotetizado que algunas arañas, como Gasteracantha cancriformis, se dispersan a grandes distancias mediante la técnica de "ballooning" (es decir, utilizando su seda para interactuar con el viento). Sin embargo, no sabemos cómo el entorno y la geografía han dado forma a su diversidad genética. Por lo tanto, probamos si los Andes contribuyeron a la diversificación de G. cancriformis actuando como una barrera absoluta o permeable para la conectividad genética entre poblaciones de esta araña en lados opuestos de la cordillera. Muestreamos miles de loci a través de la distribución de la especie e implementamos análisis de genética de poblaciones, filogenéticos y de genética del paisaje. Identificamos dos grupos genéticamente distintos estructurados por los Andes Centrales, y un tercer grupo menos estructurado en los Andes del Norte que comparte ascendencia con los dos anteriores. Esta estructura se explica en gran medida por la altitud a lo largo de los Andes, que disminuye en algunas regiones, posiblemente facilitando la dispersión y el flujo genético a través de los Andes. Nuestros hallazgos apoyan que la altitud en los Andes juega un papel importante en la estructuración de las poblaciones en América del Sur, pero la fuerza de esta barrera puede ser superada por organismos con modos de dispersión a larga distancia junto con depresiones altitudinales.

Selective concentration of iron, titanium, and zirconium substrate minerals within Gregory's diverticulum, an organ unique to derived sand dollars (Echinoidea: Scutelliformes).

Gregory's diverticulum, a digestive tract structure unique to a derived group of sand dollars (Echinoidea: Scutelliformes), is filled with sand grains obtained from the substrate the animals inhabit. The simple methods of shining a bright light through a specimen or testing response to a magnet can reveal the presence of a mineral-filled diverticulum. Heavy minerals with a specific gravity of >2.9 g/cm3 are selectively concentrated inside the organ, usually at concentrations one order of magnitude, or more, greater than found in the substrate. Analyses of diverticulum content for thirteen species from nine genera, using optical mineralogy, powder X-ray diffraction, scanning electron microscopy and energy dispersive X-ray spectroscopy, as well as micro-computed tomography shows the preference for selection of five major heavy minerals: magnetite (Fe3O4), hematite (Fe2O3), ilmenite (FeTiO3), rutile (TiO2), and zircon (ZrSiO4). Minor amounts of heavy or marginally heavy amphibole, pyroxene and garnet mineral grains may also be incorporated. In general, the animals exhibit a preference for mineral grains with a specific gravity of >4.0 g/cm3, although the choice is opportunistic and the actual mix of mineral species depends on the mineral composition of the substrate. The animals also select for grain size, with mineral grains generally in the range of 50 to 150 µm, and do not appear to alter this preference during ontogeny. A comparison of analytical methods demonstrates that X-ray attenuation measured using micro-computed tomography is a reliable non-destructive method for heavy mineral quantification when supported by associated analyses of mineral grains extracted destructively from specimens or from substrate collected together with the specimens. Commonalities in the electro-chemical surface properties of the ingested minerals suggest that such characteristics play an important role in the selection process.

Impact of persistent barrier to gene flow and catastrophic events on red algae evolutionary history along the Chilean coast.

Historical vicariance events, linked to the existence of stable physical barriers to gene flow, generate concordant genetic breaks in co-distributed species while stochastic processes (e.g., costal uplift) could cause species-specific genetic breaks as a result of local strong demographic bottlenecks or extinction. In Chile, previous studies show that the area of the 30°S-33°S could correspond to a stable barrier to gene flow that have affected the genetic structure of various algae and marine invertebrates. Here we sequenced two organellar genes (COI and rbcL) in four taxonomically accepted co-distributed red seaweeds species characterized by a low dispersal potential: Mazzaella laminarioides, M. membranacea, Asterfilopsis disciplinalis, and Ahnfeltiopsis vermicularis. Our results revealed the existence of ten strongly differentiated linages in the taxa studied. Strong genetic breaks, concordant in both space and time (divergence estimated to have occurred some 2.9-12.4 million years ago), were observed between taxa distributed across the 33°S. Conversely, in the Central/South part of the Chilean coast, the localization of the genetic breaks/sub-structure observed varied widely (36°S, 38°S, 39°S, and 40°S). These results suggest that a major historical vicariance event has modeled the genetic structure of several Chilean marine organisms in the north of the Chilean coast during the mid-Miocene, while more recent stochastic events and genetic drift could be the driving forces of genetic divergence/structuration in the central-southern part of the coast.

Genetic diversities in wild and cultivated populations of the two closely-related medical plants species, Tripterygium Wilfordii and T. Hypoglaucum (Celastraceae).

BACKGROUND: The sustainable supply of medicinal plants is important, and cultivating and domesticating them has been suggested as an optimal strategy. However, this can lead to a loss of genetic diversity. Tripterygium wilfordii Hook. f. is a medicinal plant commonly used in traditional Chinese medicine, but its wild populations are dwindling due to excessive harvesting. To protect the species and meet the increasing demand, it is urgent to cultivate it on a large scale. However, distinguishing between T. wilfordii and T. hypoglaucum, two similar species with different medicinal properties, is challenging. Therefore, it is crucial to understand the genetic diversity and population structure of these species for their sustainable utilization. RESULTS: In this study, we investigated the genetic diversity and population structure of the two traditional medicinal semiwoody vines plant species, Tripterygium wilfordii and T. hypoglaucum, including wild and cultivated populations using chloroplast DNA (cpDNA) sequences and microsatellite loci. Our results indicated that the two species maintain a high level of genetic divergence, indicating possible genetic bases for the different contents of bioactive compounds of the two species. T. wilfordii showed lower genetic diversity and less subdivided population structures of both markers than T. hypoglaucum. The potential factors in shaping these interesting differences might be differentiated pollen-to-seed migration rates, interbreeding, and history of population divergence. Analyses of cpDNA and microsatellite loci supported that the two species are genetically distinct entities. In addition, a significant reduction of genetic diversity was observed for cultivated populations of the two species, which mainly resulted from the small initial population size and propagated vegetative practice during their cultivation. CONCLUSION: Our findings indicate significant genetic divergence between T. wilfordii and T. hypoglaucum. The genetic diversity and population structure analyses provide important insights into the sustainable cultivation and utilization of these medicinal plants. Accurate identification and conservation efforts are necessary for both species to ensure the safety and effectiveness of crude drug use. Our study also highlighted the importance of combined analyses of different DNA markers in addressing population genetics of medicinal plants because of the contrasts of inheritance and rates of gene flow. Large-scale cultivation programs should consider preserving genetic diversity to enhance the long-term sustainability of T. wilfordii and T. hypoglaucum. Our study proposed that some populations showed higher genetic diversity and distinctness, which can be considered with priority for conservation and as the sources for future breeding and genetic improvement.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America.

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Phylogeographic analysis reveals extensive genetic variation of native grass Elymus nutans (Poaceae) on the Qinghai-Tibetan plateau.

Introduction: Elymus nutans holds ecological and pastoral significance due to its adaptability and nutritional value, the Qinghai-Tibet Plateau (QTP) is a key hub for its genetic diversity. To conserve and harness its genetic resources in highland ecosystems, a thorough assessment is vital. However, a comprehensive phylogeographic exploration of E. nutans is lacking. The objective of this study was to unravel the genetic diversity, adaptation, and phylogenetics of E. nutans populations. Methods: Encompassing 361 individuals across 35 populations, the species' genetic landscape and dynamic responses to diverse environments were decoded by using four chloroplast DNA (cpDNA) sequences and nine microsatellite markers derived from the transcriptome. Results and discussion: This study unveiled a notable degree of genetic diversity in E. nutans populations at nuclear (I = 0.46, He = 0.32) and plastid DNA levels (Hd = 0.805, π = 0.67). Analysis via AMOVA highlighted genetic variation predominantly within populations. Despite limited isolation by distance (IBD), the Mekong-Salween Divide (MSD) emerged as a significant factor influencing genetic differentiation and conserving diversity. Furthermore, correlations were established between external environmental factors and effective alleles of three EST-SSRs (EN5, EN57 and EN80), potentially linked to glutathione S-transferases T1 or hypothetical proteins, affecting adaptation. This study deepens the understanding of the intricate relationship between genetic diversity, adaptation, and environmental factors within E. nutans populations on the QTP. The findings shed light on the species' evolutionary responses to diverse ecological conditions and contribute to a broader comprehension of plant adaptation mechanisms.

Phylogeographic Reconstruction to Trace the Source Population of Asian Giant Hornet Caught in Nanaimo in Canada and Blaine in the USA.

The Asian giant hornet, Vespa mandarinia, is an invasive species that could potentially destroy the local honeybee industry in North America. It has been observed to nest in the coastal regions of British Columbia in Canada and Washington State in the USA. What is the source population of the immigrant hornets? The identification of the source population can shed light not only on the route of immigration but also on the similarity between the native habitat and the potential new habitat in the Pacific Northwest. We analyzed mitochondrial COX1 sequences of specimens sampled from multiple populations in China, the Republic of Korea, Japan, and the Russian Far East. V. mandarinia exhibits phylogeographic patterns, forming monophyletic clades for 16 specimens from China, six specimens from the Republic of Korea, and two specimens from Japan. The two mitochondrial COX1 sequences from Nanaimo, British Columbia, are identical to the two sequences from Japan. The COX1 sequence from Blaine, Washington State, clustered with those from the Republic of Korea and is identical to one sequence from the Republic of Korea. Our geophylogeny, which allows visualization of genetic variation over time and space, provides evolutionary insights on the evolution and speciation of three closely related vespine species (V. tropica, V. soror, and V. mandarinia), with the speciation events associated with the expansion of the distribution to the north.